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ABSTRACT 

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit 
comet-like activity by the Pan-STARRSl survey telescope using automated point-spread-function 
analyses performed by PSl's Moving Object Processing System. Deep follow-up observations show 
both a short (~ 10") antisolar dust tail and a longer (~ 60") dust trail aligned with the object's orbit 
plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and 
other well-established comets, implying the action of a long-lived, sublimation-driven emission event. 
Photometry showing the brightness of the near-nucleus coma remaining constant over ~ 30 days 
provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, 
and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we 
find an upper limit CN production rate of Qcn < 1-3 x 10^* mol s~^, from which we infer a water 
production rate of QhsO < 10^^ mol s~^. We also find an approximately linear optical spectral 
slope of 7.2%/lOOOA, similar to other cometary dust comae. Numerical simulations indicate that 
P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family 
around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m s~^. At 70 m s^^, 
this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family 
to which another main-belt comet, 133P/Elst-Pizarro, belongs. 
Subject headings: comets: general — minor planets, asteroids 
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Main-belt comets (MBCs) exhibit cometary activity 
indicative of ice sublimation yet are dy namically in- 



distinguish able from main-belt asteroids ( Hsieh fc Je 
witt||2006). Much of the current interest m studymg 



MBCs lies in the possible role of main-belt obje cts in 
the primordial de livery of terrestrial water (e.g., |Mor- 



|bidelh et al.|]2000|). Of the first five known MBCs, one 
(T76P/LINEAR) was discovered via a targeted search 
( |Hsieh ^ ^2009 j ), while the other four were discovered 



serendlpitously or by untargeted surveys (|Elst et al.|1996 
Read et al.|2"005l|Garradd et al.|2008||lNomen et al.|!2010[ ) 
Using MBCs as tracers of ice in the asteroid belt to ascer 
tain its potential for water delivery will require a much 
larger sample of these objects, and as such, dis covering 



more MBCs, ideally in untargeted su rveys (e.g., [Gilbert 



fc Wiegert||2009 Sonnett et al. 2011), is a high priority. 



Recently, main-belt objects P/^OKj A2 (LINEAR) and 
(596) Scheila have exhibited comet-like dust emission, 
but the se events a re likely due to i mpact-generated ejecta 



clouds dJewitt et al. 2010, 2011; Snodgr ass et al.||2010. 
Bodewits et"ll.|.2011. Yang fc Hsieh.,20111 [Ishiguro et| 



aT 2011 Hainaiit et 'al. 2012| . As such, these objects 



are better ch aracterized as disrupted asteroids. |Hsieh| 



et al. (2012a) considered the problem of distinguishing 
MBCs and disrupted asteroids, and concluded that re- 
current activity separated by periods of inactivity to 
be the strongest observable indicator that activity is 
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sublimation-driven, particularly when the timing of ac- 
tive episodes corresponds closely to an object' s orbital 
period (d. 133P/Elst -Pizarro and 238P/Read; |Hsieh et 
ar]|2004 2010 2011b I. Evaluating this criterion requires 
monitoring an object for many years, however. When 
observati ons of only one active episode for an object are 



available, Hsieh et al. (2012a) suggested that steady or 



increasing activity and morphological indicators could be 
used as preliminary evidence of sublimation. Direct spec- 
troscopic detection of sublimation products (i.e., gas) in 
an MBC would immediately confirm its cometary na- 
ture. Given the difficulty of obtaining sufficiently high- 
quality spectroscopy of such distant and weakly-active 
comets and of timing observations to coincide with peak 
gas p roduction though the absence of gas detections to 
date ( [Jewitt et al.|2009| [Hsieh et al.|2012"b| should not be 
considered conhrmation of the absence ot sublimation. 

2. OBSERVATIONS 

Observations of activity in P/2006 VW139 were first 
obtained by the 1.8 m Pan-STARRSl (PSl) wide-field 
synoptic survey telescope on Haleakala. PSl employs a 
3.2° X 3.2° 1.4 gigapixel camera, consisting of a mosaic of 
60 orthogonal transfer arrays, each comprising 64 590 x 
598 pixel CCDs, and Sloan Digital Sky Survey ( SDSS) 
i'- and z^-b and-like filters designated ipi and zpi (Tonry 
et al 



2012). 



ii'ollow-up imaging was performed using the 2.0 m 
Faulkes Telescope North (FTN) on Haleakala and 
Faulkes Telescope South (FTS) at Siding Spring, the Uni- 
versity of Hawaii (UH) 2.2 m telescope on Mauna Kea, 
the 1.8 m Perkins Telescope (PT) at Lowell Observa- 
tory, the 2.0 m Himalayan Chandra Telescope (HCT) at 
the Indian Astronomical Observatory on Mt. Saraswati, 
the 4.2 m William Herschel Telescope (WHT; Program 
SW2011b20) at La Palma, the 3.54 m New Technology 
Telescope (NTT; Program 185.C-1033(K)) operated by 
the European Southern Observatory (ESO) at La Silla, 
and the Lulin One-meter Telescope (LOT) in Taiwan. 
We employed 4096 x 4096 pixel Fairchild CCDs and ei- 
ther SDSS or Bessell filters for Faulkes observations, a 
2048 X 2048 pixel Textronix CCD and Kron-Cousins fil- 
ters for UH observations, the Perkins Relmaging System 
(PRISM) and Kron-Cousins filters for Perkins observa- 
tions, a 2048 X 2048 pixel E2V CCD and Bessell fihers 
for HCT obs ervations, the WH T's auxiliary-port cam- 
era (ACAM; |Benn et al .||2008| and SDSS filters on the 
WHT, the ES O l^'ainTD bjcct Spectrograph and Cam- 
era (EFOSC2; [Buzzoni et al..l98 4) and Bessell filters on 
the N TT, and a VersArray:l300B CCD ( Kinoshita et al. 
20051) and Bessell-like filters on the LOT. 



We performed standard bias subtraction and flat-field 
reduction (using dithered twilight sky images) for all 
data, except those from PSl, using Image Reduction 
and Analysis Facility (IRAF) software. PSl data were 
reduc ed using the sy stem's Image Pro cessing P i peline 



(IPP; Magnier 2006). Photometry of Landolt ( 1992 1 



standard stars and held stars was performed by mea- 
suring net fluxes within circular apertures, with back- 
ground sampled from surrounding circular annuli. For 
data obtained under non-photometric conditions, abso- 
lute calibrat ion was accomplishe d using SDSS held star 
magnitudes ( Aihara et al.|[20TT ). Conversion of r'-band 
PSl, FTN, F'i'S, and W HI' photometry to i?-band was 



acc omplis hed using transformati ons derived by |Tonr y et| 
al. (2012) and by R. Lupton ( jhttp : / /www . sdss . o rg/j) 

IJomet photometry was performed using circular aper- 
tures with varying radii depending on the nightly seeing, 
where background statistics were measured in nearby, 
but non-adjacent, regions of blank sky to avoid dust con- 
tamination from the comet. 

We also obtained longslit spectroscopy with ACAM 
(3500 A to 9400 A) on the WHT, as well as with the 
Gemini Multi-Object Sp ectrographs (GMOS; 3600 A to 



9400 A; [Hook et al.|2004[ ) on the 8-m Gemini North (GN; 
Program GN-2011B-Q-16) and Gemini South (GS; Pro- 
gram GS-2011B-Q-51) observatories. For ACAM obser- 
vations, we employed a 2'.'0-wide slit, giving a spectral 
resolution oi R ^ 300, aligned with the object's dust tail. 
For GMOS observations, we employed B600 dispersers, 
2x2 binning, and I'.'O-wide slits, giving R ^ 850, aligned 
with the dust trail. Reduction of GMOS data was per- 
formed using an IRAF package provided by Gemini. 

3. RESULTS 

3.1. Discovery of Cometary Activity 

All PSl moving object detections are screened for 
potential cometary activity by automated point-spread 
function (PSF) analyses executed nightly by PSl's Mov- 
ing Object Processing System (MOPS). This screening 
process divides the measured second PSF moment of each 
transient source by the expected PSF width as deter- 
mined from the median of all stellar PSFs in the field. 
This "psf extent" parameter is then plotted as a func- 
tion of a detection quality parameter ("psf quality") 
ranging from to 1, roughly corresponding to a normal- 
ized signal-to-noise ratio. Sources that show comet-like 
psf extent parameters and psf quality > 0.5 (to screen 
out faint sources for which measured PSF moments are 
unreliable) are flagged for human inspection and possible 
observational follow-up (Figure fl]). 

PSl observations on 2011 November 5 showed that 
the PSF of the known main-belt asteroid, (300163) 
2006 VW139 (semimajor axis, a = 3.052 AU; eccentricity, 
e = 0.201; inclination, i = 3.24°), had a FWHM o f 1"3, 
while nearby stars had PSFs with FWHMs of I'.'O (j Hsieh 

i 



et al. 2011c 



Precovery observations also obtained by 
PSl stiowed that the object exhibited a similar PSF ex- 
cess on 2011 August 30. At the time, however, the object 
was too faint (psfquality < 0.5) and was not flagged 
for inspection. In both the relatively shallow August and 
November PSl data, no extended cometary features were 
seen. Deeper follow-up observations (Table II]), however, 
revealed the presence of both an antisolar dust tail and 
an orbit-aligned dust trail (Figure[2| Section 3.2), leading 
to the object's re-designation as P72006 VW139. 

3.2. Morphological & Photometric Analysis 

The morphology observed for P/2006 VW139 from 
mid-November 2011 through January 2012 (shortly after 
perihelion) remains approximately constant (Figure ^1). 
We observe the presence of a short dust tail (~ ICr) 
pointed in the antisolar direction, a longer dust trail 
(~ 60") aligned with the object's orbit plane, and a coma 
with an intrinsic FWHM of ~ Of.'! when deconvolved with 
the seeing. This morphology is simi lar to that observe d 
for MBC P/2010 R2 (La Sagra) ( jHsieh et al.||2012b[ ). 
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as well as other w ell-established comets (e.g. 2P/Encke , 



C/Austin 1990 V; Lisse et al. 1998 Reach et al. 2000) 



Detailed dust modeling will be required to conhrm the 
nature of the observed dust features, but a simple syn- 



dyne plot (Figure [2j cf. Finson & Probstein 1968) indi 



cates that the antisolar tail is likely composed of small 
particles with short dissipation times (requiring recent 
emission to still be present) , while the orbit-aligned trail 
is composed of large particles with slow dissipation times 
(likely requiring months to traverse the observed distance 
from the nucleus). The contemporaneous observation of 
both recent and old dust emission implies the action of 
a prolonged event, consistent with sublimation and in- 
consistent with an impulsive impact. We furthermore 
compute that particles with radii of = 0.1 /xm should 
move 10" from the nucleus in just 22 days, creating an 
observable gap between the antisolar tail and the nucleus. 
No such gap is observed between November and January, 
implying continuous replenishment of the tail over this 
period, again consistent with cometary activity. 

Disrupted asteroid (596) Scheila exhibited a 30% de - 
cline in coma brightness in 8 days (Jewitt et al. 2011). 
Meanwhile, other established MBCs exhibited steady or 
incr easing coma bright ness over longer periods of time 



Hsieh et al. 2012b). Photometry of P/2006 VW 



139 



(cf. 

(Table [l]) shows no decline in near-nucleus brightness 
over 29aays between November and early December 
(though does show a 40% decline over 34 days between 
December and January). Coupled with our morpholog- 
ical analysis, the coma's photometric behavior over this 
period provides further support for the observed activ- 
ity's sublimation- driven nature. 

For reference, we compute the dust-to-nucleus scat- 
tering sur face area ratio, Ad /Aj^, determined for other 
MBCs (cf. Hsieh et al.|2012b ). Lacking independent con 
straints at this time, we use the abso lute R-hand magni- 



tude, Hn = 16.4 mag, computed by Spahr et al. (20091 
only using observations obtained well before activity was 
discovered, and find Ad/ Am ^ 2.5 during the object's 
peak observed activity in November. This level of dust 
production relative to nucleus size is comparable to that 
of 133P and 176P, but an order of magnitude lower than 
for oth er MBCs, 238P, P/2 008 Rl (Garradd), and P/La 
Sagra 
tions ( 



Hsieh et al.||2012b[ ). Multi-filter NTT observa- 
.'able [l]) indicate that the near-nucleus coma has 
an approximately Solar color of B — R = 1.06 ± 0.04, 



2009a 2010 2011a 20125) 



similar to other MBCs (Jewitt et al. 2009 Hsieh et al. 



3.3. Spectroscopic Analysis 

WHT spectroscopy was extracted using an aperture 
10" long along the slit centered on the target. A spec- 
trophotomctric standard, Feige 15, and a G2V solar ana- 
log star, HD28099, were observed to allow removal of 
atmospheric absorption features and calculation of the 
relative reflectance spectrum. The resulting spectrum 
(Figure |3^) is approximately linear with a red slope of 
7.2%/l OOOA, similar to th at seen in other cometary dust 
comae ( [Kolokolova et al.||2004[ ). 

Gemini spectroscopy was extracted using an aper- 
ture 5"1 long along the slit. A G5V solar analog star, 
HD10097, was observed for approximate flux calibration 
and to calculate the relative reflectance spectra. The 



resulting GN and GS spectra (Figure [3p) are very sim- 
ilar. Given its slightly higher signal-to-noise ratio, how- 
ever, we use the GN spectrum to compute gas production 
rates. 

We measure the standard errors in three wavelength 
regions 70 A in width and adopt the largest a of 0.057 
as the observational uncertainty of the violet_CN emis- 
sion band. We then employ a simple Haser | (l957|) 
mo del to derive the GN production rate (cf . Hsieh e 



al. 2012a Jewitt & Guilbert-Lepoutre 
resonance tluorescence effic iency ot g \ 
10~^^ erg s~^ molecule 



2012 ), using a 
ATJr= 3.63 X 
( |Schleicher| 120101 ■ We find 
an upper limit to the CN production rate of Qcn < 
1.3 X lO^'' mol s^^. No physical constraints on the CN to 
water production rate ratio in MBCs are currently avail- 
able, but adopting the average ratio in other observed 



comets (log 



Hsieh et al 



Q^n/Q ohI = -2.5; Qoh/Qh.o = 90%) (cf. 



2012a I, we infer a water production rate of 



)h20 < 10^^ mol s 

3.4. Dynamical Analysis 

To ascertain whether P/2006 VWiag is likely to be 
native to its current location, we analyze its dyna mical 
stability (cf. [Jewitt et al.|2009[ [Hsieh et al.|2012a|b[ ). We 
produce nine randomly-generated sets of 100 Gaussian- 
distributed massless test particles in orbital element 
space, centered on the object's current osculating orbital 
elements, where three sets are characterized by a val- 
ues of Ga = 0.0001 AU, CTe = 0.0001, and a, = 0.001°. 
These a values are chosen to adequately explore the sta- 
bility of the region. As a numbered asteroid, P/2006 
VWiag's true orbital element uncertainties are actually 
^ 100 — 1000 X smaller. Another three sets are char- 
acterized by a values 10 x as large, and the final three 
sets are characterized by a values 100 x as large. Treat- 
ing the eight major planets as massive particles, we then 
use the N -body integration package. Mercury (Cham- 
bers[[1999 ), to integrate the orbits of these test particles 
forward in time for 100 Myr. Limitations on comput- 
ing resources unfortunately prevent us from conducting 
significantly longer integrations (e.g., > 1 Gyr) in a rea- 
sonable amount of time. 

We find that 5% of the 1-cr test particles, 3% of the 10- 
a test particles, and 10% of the lOO-cr test particles ex- 
ceed a heliocentric distance of 50 AU, and are effectively 
ejected from the asteroid belt, within 100 Myrs (Fig- 
ureE]). P/2006 VW139 itself is found to be stable over the 
lOO^yr test period (neglecting non-gravitational forces), 
consistent with its Lyapunov time of Iiyap > 1 Myr, 
where a body is considered stable if Tiyap > 100 kyr (Tsi- 
ganis et al. 2003). No systematic distribution is evident 
for ejected l-cr or 10-cr test particles, though 44% of 100- 
cr test particles with a < 3.04 AU are ejected, perhaps 
due to the 9J:4A mean-motion resonance with Jupiter 
at 3.029 AU (Figure [IJd) . For reference, P/2006 VW139 
fluctuates between a — 3.047 AU and a — 3.069 AU over 
the course of our simulations. 
As two MBCs (133P and 176P) belong to the 



~2.5 Gyr-old (Nesvorny et al. 2003) Themis aster- 
oid family, and 133P additionally belongs to the 
young < 10 Myr Beagle sub-family (Nesvorny et al 
2008 1 , we perform 



a search for any 
tions that P/2006 VW139 may have. 



family associa- 
Employing Hi- 
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erarc hical Clustering Method (HCM; Zappala et al 
19901 analysis and analytically-determined proper or- 
bital elements retrieved on 2011 December 1 from 
AstDyS (http : / /hamilton. dm.unipi . it/astdys/), we 
compute the number of asteroids dynamically linJied to 
P/2006 VWi3g as function of cut-off distance, dc (in 
velocity space). At dc — 63 m s~^, we find a statis- 
tically significant clustering of 24 asteroids, which then 
merges with the Themis family at dc — 70 m s~^ (Fig- 
ure |4fc). This possible P/2006 VW139 sub- family is sep- 
arated from the main Themis family by several two- and 
three-body mean-motion resonances (11J:5A, 3J-2S-1A, 
and 1J-I-3S-1A) with Jupiter and Saturn, though our dy- 
namical simulations indicate these resonances may be 
only weakly destabihzing (Figure]^). P/2006 VWisg's 
association with the Themis family is interesting because 
of the family's aforementi oned association wi th 133P, 
176P, and possibly 238P ( jHaghighipour 2009), though 
more detailed study is needed to clarity the nature of 
this association, as well as the significance of the possi- 
ble sub-family. 



4. DISCUSSION 

While spectrosc opy did not reveal any evidence of gas 
emission (Section 3.3), morphological and photometric 
analyses (Section 3T5 ) strongly suggest that the activity 
of P/2006 VW139 IS cometary in nature, and not impact- 
driven. We conclude that this object is likely a MBC, 
making it the sixth such object discovered to date. The 
discovery circumstances of the known MBCs suggest that 
man y more mus t exist, particularly in the outer main 
belt f Hsieh|2009 ). Therefore, we expect that current and 
next- generation all-sky surveys like Pan-STARRS will re- 
veal more of the true extent of the population in the 
coming years, provided that techniques can be deployed 
to efficiently screen for such objects amid the enormous 
amounts of data produced by these surveys. Thorough 
physical and dynamical investigations of each new dis- 
covered object will then be essential for first determin- 
ing the most likely cause of th e observed activity (cf. 



Hsieh et al. 2012a Jewittl 2012 ), and then ascertaining 



the global properties of the population of objects con 
firmed as MBCs to better understand their implications 
for understanding terrestrial water delivery. 

Additional observations of P/2006 VW139 itself are en- 
couraged until the end (in March 2012) of the 2011- 
2012 observing window to monitor the decline of ac- 



tivity for com parison to other MBCs (cf. Hsieh et al. 
2011b| |2012b|). Observations during the 2012-2013 ob 



serving window, when the object may be largely inactive, 
should be useful for determining its nucleus properties. 
By adding to o ur knowledge of the properties of inactiv e 
MBC nuclei (cf. |Hsieh et al |2009b||Licandro et al.|2011[ ), 
the latter observations will help us better understand the 
relationship of active MBCs to the inactive main-belt 
population, perhaps facilitating the development of more 
powerful search methods, such as targeted monitoring of 
extremely likely MBC candidates. Furthermore, robust 
protocols for identifying icy objects, even in the absence 
of activity (which is transient even for known MBCs), 
will be crucial for ascertaining the true abundance and 
distribution of ice in the asteroid belt. 
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* Telescope. 

^ Number of exposures. 

" Total integration time, in s. 

Heliocentric distance, in AU. 
^ Geocentric distance, in AU. 

^ Solar phase angle (Sun-object-Earth), in degrees. 
8 True anomaly, in degrees. 

^ Position angle of the antisolar vector, in degrees East of North. 

' Position angle of the negative velocity vector, in degrees East of North. 

J Orbit plane angle, in degrees. 

^ Mean apparent magnitude in specified filter. 

' Absolute R-band magnitude (at R = A = 1 AU and a = 0°), assuming solar colors and lAU H,G phase-darkening where G = 0.15, 
where the listed uncertainty is dominated by the estimated uncertainty in G. 
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Fig. 1. — (a) Flow chart detailing PSI comet screening procedures, (b) Difference image (left) and screening plot (right) 
for sample PSI observations of 48P/Johnson where comet candidates ("selected") are plotted as red dots and "unselected" 
detections (generally consisting of inactive asteroids and false detections) are plotted as gray dots, (c) Image and screening plot 
for precovery observations of P/2006 VW139 on 2011 August 30. (d) Image and screening plot for discovery observations of 
P/2006 VW139 on 2011 November 5. 
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Fig. 2. — Orbital position plot (left) and composite images (upper right) constructed from observations detailed in Table [T] 
with the Sun (black dot) at the center, and the orbits of Mercury, Venus, Earth, Mars, P/2006 VW139, and Jupiter shown as 
black lines. Perihelion (P) and aphelion (A) are marked with crosses. Plotted positions and images correspond to observations 
from (a) 2011 August 30, (b) 2011 November 5, (c) 2011 November 12-14 (November 14 UH composite image shown), (d) 2011 
November 22 - December 4 (December 4 NTT composite image shown), (e) 2011 December 16-19 (December 19 UH composite 
image shown), and (f) 2012 January 7. Each panel is 90" x 90" with the object at the center, and North (N), East (E), the 
antisolar vector (—0), and the negative heliocentric velocity vector {—v) marked with arrows. Also shown is a syndyne plot 
(lower right) for the 2011 November 22 - December 4 period for a range of dust grain radii. 




4000 5000 6000 7000 8000 

Wavelength (A) 



9000 




3700 3800 3900 

Wavelength (A) 



4000 



4100 



Fig. 3. — (a) Relative reflectance spectrum of P/2006 VW139 from the WHT. (b) Relative reflectance spectra of P/2006 VW139 
from GN (upper panel) and GS (lower panel). Shaded regions indicate the wavelength region where the CN emission band is 
expected. Red error bars show 1-a uncertainties in the CN band region and in adjacent wavelength regions. 
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Fig. 4. — (a) Plots of semimajor axis versus eccentricity (top) and inclination (bottom) showing initial oscu lating elements of 
three 1 — a sets of Gaussian-distributed test particles whose orbits are integrated for 100 Myr (Section 3.41, with the current 
osculating orbital elements of P/2006 VW139 shown as a red circle at the center of each panel. Particles ejected in <20 Myr, 
between 20 Myr and 50 Myr, and between 50 Myr and 100 Myr are plotted with orange X symbols, yellow squares, and purple 
triangles, respectively, while particles that are not ejected within 100 Myr are marked with grey dots, (b) Same as (a) but 
for 100 — a test particle sets, where the positions of the 9J:4A, 1J-I-3S-1A, 11J:5A, and 3J-2S-1A mean motion resonances are 
marked with dotted lines, (c) Plots of proper semimajor axis vs. proper eccentricity (left) and proper inclination (right) of 
main-belt asteroids (small black dots), Themis family members (small gray dots), P/2006 VWiag (blue X symbol), and other 
MBCs (red circles). 



